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Abstract 

We propose a methodology to derive a black-hole mass for super-critical accretion flow. Here, we use 
the extended disk blackbody (extended DBB) model, a fitting model in which the effective temperature 
profile obeys the relation T e e ex r~ p , with r being the disk radius and p being treated as a fitting parameter. 
We first numerically calculate the theoretical flow structure and its spectra for a given black-hole mass, 
M, and accretion rate, M. Through fitting to the theoretical spectra by the extended DBB model, we 
can estimate the black- hole mass, M x , assuming that the innermost disk radius is r ln = 3r g ((x M x ), where 
r g is the Schwarzschild radius. We find, however, that the estimated mass deviates from that adopted in 
the spectral calculations, M, even for low-M cases. We also find that the deviations can be eliminated by 
introducing a new correction for the innermost radius. Using this correction, we calculate mass correction 
factors, M/M x , in the super-critical regimes for some sets of M and M, finding that a mass correction 
factor ranges between M/M x ~ 1.2 - 1.6. The higher is M, the larger does the mass correction factor 
tend to be. Since the correction is relatively small, we can safely conclude that the black holes in ULXs 
which Vierdayanti et al. (2006, PASJ, 58, 915) analyzed are stellar-mass black holes with the mass being 
< 100M Q . 
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1. Introduction 

The most intriguing question regarding the ultralumi- 
nous X-ray sources (ULXs) at present would be how mas- 
sive is the black hole in each ULX. 

Known to be very bright extragalactic X-ray sources 
[L x ~ 10 39_41 ergs _1 ), brighter than the ever-known 
Galactic X-ray sources, ULXs are good candidate to prove 
the existence of the hypothetical new class of black holes, 
the intermediate mass black holes (IMBHs) , whose masses 
range over 10 2 - 10 4 M Q (e.g. Colbert & Mushotzky 1999; 
Makishima et al. 2000). If we simply assume L x < Le, 
where Le is the Eddington luminosity, M > 10 3 M Q for 
L x ~ 10 41 erg s _1 . This provides strong support to the 
IMBH hypothesis. 

The success in fitting the spectra of several ULXs with 
a multicolor disk blackbody (DBB) and a power-law (PL) 
model further supports the idea of IMBHs for the ULXs, 
since the cool inner disk temperatures and the large inner- 
most radii, obtained through model fitting, suggest a black 
hole mass within the IMBH range (Miller et al. 2003,2004; 
Cropper et al. 2004, Roberts et al. 2005). Moreover, the 
explanation behind the IMBH notion for ULXs is based 
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on the standard accretion disk model (Shakura & Sunyaev 
1973), a well-established model for high-energy release 
mechanism through gas accretion flows around a compact 
object. 

On the other hand, some attempts have been made 
to challenge the IMBH interpretation of ULXs (King et 
al. 2001). Several models have been proposed to allow 
stellar-mass black holes to power ULXs. In fact, such at- 
tempts are very reasonable since the existence of stellar- 
mass black holes is well founded while IMBH formation 
has not been well understood (e.g. Madhusudhan et al. 
2006). 

Amongst several models proposed, a model of super- 
critical accretion flow onto a stellar-mass black hole has 
been quite favorable recently (e.g., Watarai et al. 2001; 
Okajima et al. 2006; Vierdayanti et al. 2006; Tsunoda 
et al. 2006). One attractive model supporting a super- 
critical accretion flow onto a stellar-mass black hole is the 
slim disk model (Watarai et al. 2000; Kawaguchi 2003, 
Ebisawa et al. 2003; Okajima et al. 2006) introduced 
in the late 1980s by the Warsaw and the Kyoto groups 
(Abramowicz et al. 1988, 1989, see Kato et al. 2008 for a 
review). Its application to some Galactic black holes was 
investigated by Watarai et al. (2000) in their attempts to 
explain the nature of some bright X-ray source binaries. 
Later on, Watarai et al. (2001) also proposed the slim 
disk model for ULXs. 
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The black-hole mass mentioned above was estimated 
from an X-ray spectral fitting in which the innermost ra- 
dius, ri n , was assumed to coincide with that of the last sta- 
ble circular orbit. r ms (= 3r g for a non-rotating black hole, 
where r g [= 2GM/c 2 ] is the Schwarzschild radius). This 
assumption has been proven to be good for the cases of 
Galactic black hole binaries (e.g., Makishima et al. 1986; 
Tanaka & Shibazaki 1996; Ebisawa 1999; Dotani et al. 
1997). 

However, numerical calculations have shown that as 
mass accretion rate increases, the flow inside 3r g might 
produce significant emission (Watarai et al. 2000). The 
heat generated near 3r g is trapped by the inward mate- 
rial flows, and is emitted inside 3r g . Some amount of the 
trapped radiation going closer to the black hole is swal- 
lowed. Therefore, since the situation is more complicated 
in high mass accretion rate systems, we cannot simply re- 
late the innermost radius from the fitting with the last sta- 
ble circular orbit. Indeed, Watarai and Mineshige (2003a) 
have stressed that the observed value of the inner edge of 
the disk is more like a radiation edge, outside of which 
substantial emission is produced. In this work, we in- 
vestigated the assumption r- m = 3r g for the super-critical 
accretion flows around a non-rotating black hole and cal- 
culate corrections to make more accurate black-hole mass 
estimations. 

The plan of this paper is as follows: 

We describe the disk structure and spectrum calcula- 
tions in the first two subsections in section 2, and the 
fitting model, together with the conventional method to 
derive black hole mass from the fitting in the remaining 
subsection. 

In section 3, we present our fitting analysis and results, 
while focusing on some sets of parameters. 

Section 4 is devoted to a discussion, and finally section 
5 concludes the paper. 

2. Methodology 

To avoid confusion, we will first describe the overall 
picture of our methodology which is also summarized in 
table 1. The detail steps in the fitting are explained in 
another section (section 3). 

The main steps of our methodology are as follows: We 
first calculate the disk structure for various values of 
black-hole mass, M, and the mass-accretion rate, M (sub- 
section 2.1). Next, we use the outputs from the disk struc- 
ture calculation to obtain the spectra (subsection 2.2). We 
follow the calculation by Watarai et al. (2005) for both 
steps. We then fit those spectra with the extended DBB 
model and derive the black hole mass, M x from value of 
r m obtained from the fitting (subsection 2.3). That is, we 
assume r- m = 3r g as that of the DBB model. Finally, we 
obtain the correction factor, M/M x , the ratio of the actual 
black hole mass to the derived black hole mass (subsection 
3.2). 



2.1. Disk Structure Calculation 

The standard accretion-disk model has been very suc- 
cessful in describing optically thick flow structure as long 
as the mass accretion rate of the systems in question is less 
than the critical mass-accretion rate 1 , M < M C rit = Le/c 2 . 

It is possible, however, that the mass-accretion rate of 
a system exceeds the critical mass-accretion rate in disk 
accretion because of anisotropic radiation field (Ohsuga 
et al. 2005). A disk in which its accretion rate exceeds 
the critical value is called a super-critical accretion disk. 

In this section we describe the calculation of the struc- 
ture of such a disk (e.g., Watarai et al. 2000, 2005; see 
reviews by Kato et al. 2008). The equations are written 
by using cylindrical coordinates (r, ip, z). A non-rotating 
black hole is assumed as the central object and a pseudo- 
Newtonian potential (Paczyhski & Wiita 1980) is adopted 
to model the general relativistic effects in the gravitational 
field of the black hole, 



GM 

R-r e 



(1) 



where R = \/r 2 + z 2 . 

The pressure and density are related to each other by 
polytropic relation (Hoshi 1977), poa p 1+ ~ , in the vertical 
direction, where N is the polytropic index. The vertically 
integrated density, E, and pressure, II, are then 
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and H is the scale height. 

Assuming hydrostatic balance in the vertical direction 
we have 

Po 
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where Ok is the Keplerian angular velocity of the disk ro- 
tation in the pseudo-Newtonian potential, which is defined 

by 
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Note that in terms of II and S, the hydrostatic balance in 
the vertical direction can be written as 



Qk£T 2 = (27V + 3) 



n 



(7) 



Regarding the polytropic index, we set N = 3 for the whole 
disk calculation. 



Some papers use M cr ; t = Lg/ (r]c 2 ) with rj, the energy conversion 
efficiency, being a constant. However, this may cause confusion 
since rj is not constant in the slim disk model. 
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Table 1. The overall picture of our methodology. 



Step 


Calculation 


Input 


Output 


Section 


f 


Disk structure 


M, M 


Teff(r), H(r), v r (r), v v (r) 


2.1 


2 


Spectrum 


T eff (r), H(r), v r (r), v v {r) 


Sy 


2.2 


3 


Fitting 


s v 




2.3 


4 


Correction factor 


Tin, M 


M/M x 


3.2 



The other basic equations are obtained by integration 
in the vertical direction under the assumption that the 
radial velocity, v r , and the specific angular momentum, £ 
(= rv v = r 2 fl), do not depend on the vertical coordinate. 

The continuity equation gives 

—2irrY>v r = M = constant, (8) 

where M is the mass accretion rate. 

The radial component of the momentum equation can 
be written as 



\dXi_£ 2 -£\ ndlnf2 K 



dr £ dr 



dr 
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where £k is the Keplcrian angular momentum, defined by 
£ K =r 2 Q K . The last term of the above equation, § rfl '^ 2K , 
is a correction term resulting from the fact that the radial 
component of the gravitational force changes with height 
(Matsumoto et al. 1984). 

The angular-momentum balance is obtained from the 
(^-component of momentum equation, and can be written 
as 



M(£-£ in ) = -2irr 2 T r , 



(10) 



where £- m is an integration constant, which represents the 
specific angular momentum finally swallowed by the black 
hole, and T rtp is the vertical integration of t rtpi that is 

(11) 
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t rip dz = —all, 



with t rip = —ap as in the standard Shakura-Sunyaev pre- 
scription. 

The energy equation can be written as 
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where Q„, ]v is the advective cooling, defined by 
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where s is the specific entropy on the equatorial plane. 
The viscous energy generation rate is given by 
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while the radiative cooling rate is given by 

ad dz = 2F, (15) 



Qrad = / Irs 
J-H 



where F is the radiative flux per unit surface area on the 
disk surface, given by 
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where a is radiative constant, T c is the temperature on 
the equatorial plane, and r is given by 

(17) 
(18) 

(19) 



r = kE = (k cs + Kff)S. 
The equation of state is given by 

n = n gas + n ra d = _ B -^t,t c + n ra d, 

where II ra( j is defined by 

f H 1 / z 2 \ 4 

^Lr^v-jp) dz - 



We first specify the values of M, a, M, and an initial 
guess of £j n . The outer boundary conditions arc imposed 
at r = 1 X 10 4 r g , where all physical quantities are taken as 
those of the standard disk. The inner calculation bound- 
ary is taken at r ~ 2.0r g . Note that it is important to solve 
the equations, even inside the marginally stable last cir- 
cular orbit at 3r g , since substantial amount of matter still 
exists there. In fact, the emission from this region makes 
a rather important contribution when M > 100A/ cr it. Also 
note that it is possible to calculate the disk structure in- 
side 2r"g, but in this region the disk in fact becomes op- 
tically thin and thus the blackbody assumption made in 
calculating the spectrum may no longer be valid. Further, 
the contribution from the innermost part at r < 2r g is 
small because of a photon redshift. Thus, we left this 
region for our future studies. 

The solution should satisfy the regularity condition at 
the transonic radius located between (2 - 3) r g and the 
free boundary conditions at the inner edge. The basic 
equations are first radially integrated for an a priori given 
£- m . When l\ a is larger (or smaller) than the correct value, 
the radial velocity vanishes (diverges) before r g . Iteration 
of is therefore needed to find an appropriate value for 
which a transonic solution is obtained. In other words, 
the value of £- ln is determined uniquely by the transonic 
condition. 

2.2. Spectrum Calculation 

In calculating the spectrum, we apply a calculation 
method proposed by Watarai et al. (2005). They 
adopted the 'Ray- Tracing method' (Luminet 1979; Fukue 
k, Yokoyama 1988) in calculating the spectra of a disk 
around a Schwarzschild black hole. Before reaching the 
observer, a photon emitted from some point on the disk 
travels along the null geodesies. In the Schwarzschild 
space-time, the photon's trajectory is determined by 



K. Vicrdayanti, K.Watarai, and S. Mincshige 



[Vol. 



dx 2 



3r g 
2T 2 ' 



(20) 



where r is the distance from the center (black hole), and 
X is the angle measured from the direction of the observer 
(see reviews by Kato et al. 2008). The trajectories will 
then satisfy the following differential equation: 
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where b is the impact parameter. 

The light ray is calculated from an observer's coordi- 
nate to the surface of the disk by the 'Runge-Kutta-Gill 
Method' and the scale-height is defined as the disk sur- 
face of the last scattering. The numerical integration is 
stopped when the ray arrives at the surface of the disk, 
and the final arrival point is interpolated from the cal- 
culated data of the disk model. The obtained quanti- 
ties, such as temperature, velocity field, and redshift on 
the disk surface, are transformed to the observer's values, 
which are used to calculate the observed spectrum. 

The observed spectra were calculated by using the 
Lorcntz invariant relation. 
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where I Vtm and z/ om are the intensity and frequency in 
the disk co- moving frame, respectively. Some simplifica- 
tions are also made here, that the blackbody radiation, 
B v (T e f{), is assumed in order to calculate the spectra, and 
the temperature on the disk surface is adopted as the effec- 
tive temperature. The thermal Comptonization process, 
which may become important in the innermost region, is 
approximated by introducing a spectral hardening factor, 
K. Wc set k= 1.7 throughout the calculation (see discus- 
sion in subsection 4.2). The local flux is calculated by 
using a diluted blackbody approximation, 

^jTtB u {kT cS ). (23) 
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The redshift factor, (1 + z), which consists of a gravita- 
tional redshift part and the Dopplcr parts, is given by, 
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where L, 7, and D are the lapse functions representing the 
gravitational redshift, Lorentz factor, and Doppler factor, 
respectively (Kato et al. 2008). 

According to Watarai et al. (2005), the inclination- 
angle dependence of the accretion disk spectrum is com- 
posed of the following three factors: 

1. The projection effect, which is the change of the ef- 
fective area with the change of the inclination angle. 

2. The effect of Doppler beaming towards the ob- 
server's direction, which causes asymmetry in the 
flux distribution with respect to the rotation axis. 

3. The self-occultation effect (Fukuc 2000), that is, the 
central part of the disk is obscured by the cooler 



outer part of the disk. This effect is significant for 
high-inclination systems. 

Note that self-irradiation (Cunningham 1975) was not 
considered in the present study. 

The spectrum calculation done by Watarai et al. (2005) 
was not intended for calculating the spectrum over a large 
area of the disk. We made a simple modification to make 
the calculation for larger disk area possible with reason- 
able calculation time. This is necessary since the spectra 
of the high mass-accretion rate disks extend into the soft 
X-ray region (Szuszkiewicz et al. 1996). To do so, we di- 
vided the disk by following y/r pattern. By doing so, we 
can keep the detail of the innermost part without losing 
those of the outermost part too much. 

2.3. Fitting Model 

In this subsection we switch to a description of the fit- 
ting model and the conventional method to derive the 
black hole mass. 

2.3.1. Extended disk blackbody model 

The fitting model that we describe in this section, and 
is the main tool used in this work, is the extended disk 
blackbody (extended DBB) model. This X-ray spectral 
fitting model was originally proposed by Mincshige et al. 
(1994) to investigate whether the temperature gradient 
found from the observations of Nova Muscae 1991 could 
deviate from the canonical value of the standard disk. 0.75 
(T cff cxr- 3 / 4 ). 

The name 'extended DBB' is chosen, because it is a 
natural extension of the DBB model 2 . It becomes a gen- 
eralization and an extension of the DBB model, since the 
temperature gradient is now treated as one of the fitting 
parameters. The temperature profile, hence, becomes, 



T. 



off : 
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where T- ln and n n have the same interpretation as those 
in the DBB model. 

By using self-similar solutions, Watarai and Fukuc 
(1999) found that the temperature distribution of the 
slim disk is somewhat flatter than that of the standard 
disk, that is T cff cx r" 1 / 2 (see also Wang & Zhou 1999). 
Therefore, the temperature profile of the standard disk 
is proportional to r~ 3 / 4 , while for the case of slim disk, 
r -i/2_ 

Generally speaking, hence, by using the extended DBB 
model, r~ p , for the fitting, we can distinguish the case 
where the standard disk is more plausible than the slim 
disk, and vice versa, from the spectrum profile. 
2.3.2. Conventional methodology to derive black-hole 
mass 

When we obtain a good fit to the observed spectra by 
the DBB model, we can easily estimate the black-hole 
mass and the Eddington ratio, L/Le (with L and Le be- 
ing disk luminosity and the Eddington luminosity, respec- 
tively), based on the standard disk theory 



In some papers it is called the 'p-frce'model, but we do not use 
it, since this is a confusing terminology. 
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Fig. 1. Comparison between the standard disk model (Shakura and Sunyaev 1973; analytical Newtonian non-transonic solution) in 
solid lines and the numerical model (pseudo-Newtonian, transonic solution) at low m (m = 1), in dashed lines. Left panel: effective 
temperature profiles. Right panel: Spectra. We also plot analytical pseudo-Newtonian non-transonic solution in thin dot-dashed 
lines. 



The basic methodology is summarized in Makishima et 
al. (2000). 

We write the bolometric luminosity of an optically thick 
accretion disk as 

ibo^^jW^)- 1 , (26) 

with /boi being the bolometric flux, i being the inclination 
of the disk (i = corresponds to face-on geometry) and D 
is the distance. This Lboi is related to the maximum disk 
color temperature, Tj n , and the innermost disk radius, rj n , 
as 

£boi = 4 7 r(r in /0VT in // t ) 4 , (27) 

where k ~ 1.7 (Shimura & Takahara 1995) is the ratio 
of the color temperature to the effective temperature, or 
spectral hardening factor, and £ = 0.412 is correction fac- 
tor reflecting the fact that T- ln occurs at somewhat larger 
than r- m . 

Hence, the innermost disk radius is 

We may identify n n with the radius of the last stable 
Kcplerian orbit. Thus, we may in general write 

M 1 
r- m = 3/?r g = 8 - 86 / ? (^) km with g < P < 1 (29) 

by which wc can determine the black hole mass. Note 
that (3 depends on the black hole spin and j3 = 1/6 for 
an extremely rotating Kerr black hole, and (3=1 for a 
Schwarzschild black hole. In the present study, we focus 
on the cases of the Schwarzschild black hole (/? = 1). 
2.3.3. Problem with the Conventional Methodology 

The extended DBB model is a phenomenological model 
which was proposed to investigate a possible deviation of 
the temperature gradient from 0.75, the value predicted 
by the theory of the standard disk. It turned out that 



this model can also be useful to discriminate the slim disks 
from the standard disks, as explained in the previous part. 

However, physical interpretation of the fitting results 
by using the extended DBB model is far from being clear, 
unlike in the case of DBB model. We need a methodology 
to make a physical interpretation of the extended DBB 
model in the light of the super-critical accretion-disk the- 
ory. 

In the case of ULXs, Vierdayanti et al. (2006), for ex- 
ample, have tried to make a physical interpretation of the 
fitting results from the extended DBB fitting model. They 
adopted a method that is widely used for the DBB fitting 
model which corresponds to the standard disk theory to 
derive black hole masses from the fitting results. The fit- 
ting shows the tendency that the observational data of the 
analyzed ULXs are more likely to be fitted with the slim 
disk model. We realized that the adoption of the conven- 
tional method used for the DBB model is not fully appro- 
priate, since the spectral shapes of the slim disk is different 
from those of the standard disk. Therefore, the adoption 
of the conventional method is inadequate for estimating 
the black hole mass. In addition, a number of relativis- 
tic effects as well as self-occultation effect certainly affect 
the spectra, and thus the fitting results. The implication 
of such relativistic cases are not obvious, and have been 
poorly investigated so far. In conclusion, some correction 
should be introduced to permit this adoption. To be pre- 
cise, in the following section we will try to find a mass 
correction factor, defined as m/m x , where m (= M/Mq) 
and m x (= M X /M Q ) are the original black-hole mass and 
the black-hole mass obtained from the fitting, respectively. 
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Table 2. The fitting steps (see text for definition of £ trans , and M/M x ). 



Step 


Fitted Model 


m-value 


Corrections 


Section 





Shakura-Sunyaev model 


1,10 


£ = 0.412 


3.1 


1 


Numerical transonic model 


1,10,32,100,320,1000 


£trans = Q 353^ M/M* ( for high-M) 


3.1 and 3.2 



c 



> 



3.5 
3 

2.5 '^b 

2 
1.5 



1.7; 

1.5: 

1.3 

1.1 

0.9 

0.58 



0.54 

0.5^ 
0.46 
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m=100 
m=320 
, m=1000 



10 



20 



40 



log m 



80 



fcT in (middle 
320, and 1000 



Fig. 2. Fitting results: p (bottom panel), 
panel), and f; n /r g (top panel) with rh = 100, 
and i = for various m. 

3. Fitting Analysis and Results 



We fit the theoretical spectra of the slim disk (Watarai 
et al. 2005) with the extended DBB model. In the the- 
oretical spectrum calculation, there are three parameters 
that should be determined in the disk structure calcu- 
lation, whose results are used to calculate the spectra. 
They are, the black-hole mass, m(= M/Mq), the mass- 
accretion rate, m(= Mc 2 /£e), and the viscosity parame- 
ter, a. Regarding the viscosity parameter, in this study 
we only focus on a = 0.01. Another parameter is needed in 
the spectrum calculation, that is the inclination angle of 
the disk or the viewing angle of the observer. On the other 
hand, there are three fitting parameters in the extended 
DBB model, T- m , r- ln , and p. We divide our fitting proce- 
dure into two important steps, which are summarized in 
table 2. 

To relate r- m with the radius of the inner edge of the 
disk, we need a correction, £. However, we found numeri- 
cally that the £-value does not coincide with the analyti- 
cally expected one, £ = 0.412. Hence, we firstly fit two disk 
models, analytical standard disk model and numerical disk 
model for low ro-case, with the extended DBB model to 
find the basic corrections, which will be explained in the 
following subsection. The details are given in subsection 



3.1 

Secondly, after we found the basic corrections we ap- 
plied them to the high m-casc. This step is described in 
subsection 3.2. 

3.1. Basic Corrections 

As a first step of our fitting, we elucidate some basic cor- 
rections that were set aside in some previous studies. We 
call them basic corrections, since they emerged from the 
disk structure calculation before we calculated the spec- 
tra. 

In the fitting with the extended DBB model, which 
is just an extension of the DBB model, we expect the 
same results for both analytical standard disk model and 
numerical disk model for the case of moderately low rh, 
rh < 1. This is because theoretically, the disk should be- 
have like a standard disk when rh is moderately low. On 
the other hand, it should deviate from that of the standard 
disk when rh is moderately high, since the assumptions 
made for the standard disk case will break down. 

Interestingly, our fitting results show some discrepan- 
cies even for a moderately small rh case. We studied some 
previous works and come to a conclusion that a possi- 
ble source of discrepancies is q small deviation of the nu- 
merical disk model from that of the standard disk. (cf. 
Ebisawa et al. 1991; Kubota et al. 1998: Kubota 2001). 
In order to prove it, we compared the numerical tran- 
sonic flow model in a pseudo-Newtonian potential with 
the analytical Shakura and Sunyaev (SS) model (in which 
the transonic nature of the flow is not considered) in a 
Newtonian potential. 

The left panel of figure 1 shows the effective temper- 
ature profiles of both models. The solid lines represent 
the analytical result of Shakura-Sunyaev model, while the 
dashed lines represent that of the transonic model for 
m — 10 and rh = 1. We can see that even at low rh, 
a discrepancy exists, especially in the innermost region, 
where the transonic solution plays its role and the effect 
of different potential forms is more significant. We then 
calculate the spectra of both models by using the spec- 
trum calculation code of Watarai et al. (2005). However, 
here and only when it is stated, we neglect the relativistic 
Dopplcr and light bending effects for simplicity. The re- 
sults are shown in the right panel of figure 1. In addition, 
we also plot analytical pseudo-Newtonian non-transonic 
solution as a comparison (thin dot-dashed lines in figure 
1). As shown in figure 1, the effect of the difference in the 
choice of potential form is more prominent compared to 
the effect of transonic/non-transonic solution. 

The discrepancy in the temperature profiles shows that 
we need a correction factor in addition to n for the in- 
ner temperature. In order to find the required correction 
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Table 3. Fitting results with the extended DBB model for 
a = 0.01 and i = 0. The unit for fcT; n is keV and energy range 
for fitting is 0.3 — 10 keV. The luminosity is calculated for all 
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factor, we calculate the spectrum of those two disk struc- 
ture models without considering any spectral hardening 
effect, i.e., we set K = 1 and compare the values of Tj n . 
We found that T m of the numerical model is by a factor of 
~ 1.2 lower than that of the standard disk. Thus, in fitting 
the numerical disk model with the extended DBB fitting 
model, the inner temperature correction /Ti n = 1.2 is 
required in addition to k, where Tj^ in the inner tem- 
perature of the Shakura-Sunyaev model. Therefore, it is 
necessary to define K rcal = K trans K, as a real correction for 
the inner temperature, where K trans = 1/1.2 and k ~ 1.7 
for Compton scattering. 

We found that the innermost radius, ?' in , also deviates 
from that of the standard disk. Note that in calculating 
ri n , we have used the bolometric luminosity, and thus the 
luminosity does not depend on the fitting energy range. 
We decided to handle this deviation separately from that 
of the temperature. The reason is to keep the black-hole 
mass estimation method simple. That is, we can simply 
obtain r; n from the Tj n (without correcting it first), and 
just add the correction at the end of the procedure, that 
is after we obtain the estimated black hole mass. The 
correction for T- ln becomes important when we want to 
estimate the mass accretion rate which is not the focus of 
our present study. 

By using the uncorrected T ln in the fitting, the deviation 
in the innermost radius can be eliminated by changing the 
correction factor, £, introduced earlier (see subsubsection 
2.3.2). In our present study, we found that the new value 
°f £> £ trans = 0.353 is more appropriate instead of £ = 0.412 
which is commonly used (cf. Gicrlihski et al. 1999). We 
summarize the fitting results in table 3. 

We notice that the p- value obtained by the fitting is not 
close to p = 0.75 but significantly lower, p < 0.64. This 
reflects a bit flatter temperature profile of the transonic 
solution (see figure 1) 

3.2. Fitting Results for High m-case 

Now, for applying to the ULXs, we examine the high in- 
case. Here, we apply the correction that we found in low 
m-case to the high m-case. That is, we use £ trans = 0.353. 
Note that in addition to K we also found necessary correc- 
tion for the inner temperature, which is K trans = 1/1.2, 
and thus the real correction for the inner temperature 
is K rcal = k/ 1.2. However, as we decided to treat this 
correction separately, the most important basic correc- 
tion for our fitting methodology is the new value of £, 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

log Energy [keV] 

Fig. 3. Fitting range sample 

which is £ trans = 0.353, while still using k ~ 1.7 instead of 
K rcai _ K trans K f Qr ^ e fitting. The detailed results of the 
fitting are summarized in table 4 and 5. 
3.2.1. m-dependence 

Firstly, let us see the fitting results for different black 
hole masses. We calculated the spectra of the slim disk 
for m = 10, 30, and 100, fixed m at 100, and considered 
the face-on case only, i — in this subsection. Let us focus 
on the inner temperature, kT- ln , and p (figure 2). 

We see that the inner temperature, fcX| n , decreases as 
the mass increases. This is natural, since the more massive 
is the black hole the lower is the temperature, as T m oc 
M -1 / 4 . Interestingly, the value of p behaves in the same 
way, which is not obvious. It is important to note that 
p is determined by the slope in the lower energy (soft 
X-ray) region (Watarai et al. 2000), that is the middle 
part of the spectrum. The disk spectrum is shifted to the 
lower energy region as the mass increases (see Fig. 3). In 
other words, a fixed fitting energy range corresponds to 
a different part of the spectrum for different black hole 
mass (see figure 3 for illustration) . This is the reason why 
p depends on the black-hole mass. 

We found that the higher is the black hole mass, the 
smaller is the p- value (see table 4). Note that recent X- 
ray fitting of some ULX spectra with the extended DBB 
model also showed that p ~ 0.5. This is natural if we 
jump to conclusion that those ULXs have high black-hole 
masses based on our finding above. However, a small p 
can also be obtained even by a stellar-mass black hole 
with a high mass-accretion rate, which is discussed more 
in the next subsection. The important thing is that the 
dependence of p on black-hole mass that we found here is 
natural, since we restrict the fitting energy range. 

Now, let us see the behavior of the inner radius obtained 
from the fitting. We found that the inner radius also de- 
pends on the black-hole mass; a larger black-hole mass 
has a smaller r- ln /r g . In the fitting, the innermost radius 
is mainly determined by the high-energy part of the spec- 
trum (inner region of the disk). The dependence of the 
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inner radius on the black-hole mass shows that the shapes 
of the spectrum from the inner part of the disk change 
with the black-hole mass. In other words, the spectrum 
of the inner part of the disk of a massive black hole is not 
a simple scale up of that of the stellar-mass black hole 
(cf. Abramowicz et al. 1988). However, the difference is 
small and becomes less prominent in the lower energy part 
(outer part of the disk). 

3.2.2. m- dependence 

Now, let us see the fitting results for various mass ac- 
cretion rate, rh (table 4). The most interesting result is 
the luminosity. In the theory of a slim disk, the luminos- 
ity at first increases as the mass-accretion rate increases. 
However, above a certain value, the radiation efficiency 
decreases as the advection effect becomes dominant, and 
the luminosity increases not at the same rate as the mass- 
accretion rate (Jaroszyhski et al. 1980). In other word, 
the luminosity is saturated above a certain mass-accretion 
rate value. 

We found that the luminosity rapidly increases as we 
increase rh from 32 to 100. As we further increase rh up 
to 1000, slower increase in luminosity is realized (figure 
4). In other words, our results are in good agreement 
with the theory, and thus justifies our methodology Note 
that the fitting results of the luminosity for different black 
hole masses coincide in figure 4. 

As we briefly mentioned above, p decreases as rh in- 
creases which is also in good agreement with theoretical 
prediction. Since as rh increases the heat is trapped in the 
flow of matter and thus the advective, horizontal heat flux 
becomes important (Abramowicz et al. 1988). In other 
words, the effective temperature gradient becomes flatter 
as rh increases, and thus p decreases. Although higher 
black hole mass gives lower p-value, it is important to re- 
member that a higher black-hole mass also gives a lower 
Tj n . Therefore, it is more likely that a lower mass black 
hole with higher mass accretion rate can explain high tem- 
perature, low p-value ULXs. 

3.2.3. i-dependence 

Lastly, we would like to show the behavior of the spectra 
for various inclination angles. We mentioned earlier when 
we described the spectrum model of Watarai et al. (2005) 
that the spectral shape of their model strongly depends on 
the inclination angle as the mass accretion rate increases. 
Let us now see the fitting results for various inclination 
angles. Here, we fixed rh = 100, and m = 10, 30, and 
100 (table 5). The inner temperature of the disk at first 
tends to increase as the inclination angle increases, but a 
sudden drop is realized at an inclination angle above 50° 
(figure 5). The slight increase in kT- m at low inclination 
angles is basically caused by the Doppler boosting, while 
a sudden drop above 50° is caused by the self-obscuration 
effect. Note that we obtained very similar trend even for 
different black hole masses. In this study, we focus on 
the case for inclination angle below 50° to explain ULXs 
with high inner temperature, as suggested by Watarai et 
al. (2005). 
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Fig. 5. Fitting results: kT in with m = 10, 30, and 100 and 
rh = 100 for various i. 

3.2.4- Mass Correction Factor, m/m K 

Let us now examine the mass correction factor, which 
wc define as the ratio between the real black hole mass, m, 
to the estimated mass by the X-ray spectral fitting, m x . 
Here, we adopt the assumption used in the conventional 
disk blackbody model, that is n n = 3r g . 

Firstly, let us see the dependence of the correction fac- 
tor on rh (figure 4). Here, we fix i = and consider three 
different masses: m — 10, 30, and 100. Interestingly, we 
can see a similar trend in all cases; the correction factor in- 
creases as the mass-accretion rate increases. Interestingly, 
we can see that the correction factor also seems to depend 
on the black-hole mass (see table 4). It follows a similar 
trend as that of rh; that is, a larger correction factor is 
needed for a larger black- hole mass, although it is less 
prominent for the small rh case. 

These findings show that correction factor (~ 2) is 
needed when we work with the system of high m and rh. 
Unfortunately, those are indeed two pieces of information 
that we intend to find. However, since we are interested 
in the high temperature ULXs, we suggest the correction 
factor of 1.2 - 1.3 and 1.3 - 1.6 for m = 10 and 30, respec- 
tively, depending on the mass-accretion rate. 

Next, let us see the dependence of the correction factor 
on the inclination angle (figure 6). Wc found that as long 
as we consider inclination angles below 50° , the correction 
factor does not have a strong dependence on the inclina- 
tion angle. Here, we show three examples for rh = 100: 
m = 10, 30, and m = 100. 

To summarize, for high m-case, we need a mass correc- 
tion factor defined as m/m x . 

Before leaving this section, it is important to note that 
in calculating the disk structure, and thus the spectrum, 
we set the inner calculation boundary at 2r g . On the 
other hand, to derived the black hole mass we adopt 
r- m = 3r g . Since the correction factor that we found here 
varies with m and also rh, we conclude that the inner ra- 
dius obtained from the fitting does not represent the inner 
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Fig. 4. Fitting results: L/Le (left panel) with j: 
calculated for m = 10, 30, and 100 with i = 0. 
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Table 4. Fitting results with the extended DBB model for 
a = 0.01 and i = 0. £ trans = 0.353 is assumed to derive r ln and 
the real fcT; n is the value in the table multiplied by re trans . 



Fig. 6. Correction factor (= m/m x ). Top panel: for 
m = 10, 30, and 100 and m = 100 for various i when i is 
unknown. Bottom panel: same as top panel but when i is 
known. 

edge of the disk as was stressed in Watarai and Mincshige 
(2003a). They suggest that the inner radius from the 
fitting represents the radiation edge outside which sub- 
stantial emission is produced. In fact, for super-critical 
accretion disks, the emission produced inside 3r g is really 
significant (Watarai et al. 2000). 

3.3. X-ray HR Diagram 

Finally, we plot our results in the so-called X-ray HR 
diagram, which shows the relation of X-ray luminosity, L x , 
with the inner temperature (in keV), fcTm, obtained from 
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the fitting (figure 7). Here, we now include the correction 
for T in in the X-ray HR diagram. That is T in = nT cS /l.2 

We plot for m= 10 (in black '+'), 30 (in black 'x '), and 
100 (in black '*') and m = 1, 10, 32, 100, 320, and 1000, 
while i = 0. This figure is similar to that of Watarai et 
al. (2005). However, in their paper, they only calculated 
for m = 10, and then interpolated the results for higher 
black-hole masses. We also have included the required 
corrections discussed in subsection 3.1. 

For low rii, our results follow L oc T^ n , as expected from 
the sub-critical accretion (standard accretion disk) the- 
ory. For high ?ti our results follow L oc T? n . We also plot 
some ULXs and well-known BHCs by using the references 
in Vierdayanti et al. (2006). For BHC whose inclina- 
tion angle has been known, we included it in calculating 
the luminosity, while for those whose inclination angle is 
unknown, we assumed a face-on case. Some ULXs can 
very likely be explained as stellar-mass black holes with a 
super-critical mass accretion rate. Except for the source 
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Table 5. Same as table 4 but here we fixed m = 100 while 
varying the inclination angle, i. 
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IC 342 SI, Vierdayanti et al. (2006) fitted the other ULX 
data with the extended DBB model. The data of IC 342 
SI are taken from Mizuno et al. (2001), in which they 
fitted the ASCA data with the multicolor disk blackbody 
(DBB) model. 

The red dotted lines shows the masses and rj (= L/Le) 
which were calculated from the standard disk relations 
[see Makishima et al. 2000 equation (9) and (11)]. We 
also plot the fitting results from Shakura and Sunyaev 
analytical solutions (see section 3.1) in red 'Cgi'mark. We 
can thus safely conclude that the black holes in ULXs, 
which Vierdayanti et al. (2006) analyzed, are stellar-mass 
black holes with mass < 1OOM . 

4. Discussion 

In this section we would like to discuss some important 
issues that arc related to our present study. The first 
and the most important issue is the dependence of the 
fitting results on the fitting energy range. The next two 
issues are the spectral hardening factor, which remains an 
open question for the case of super-critical accretion, and 
outflow of material from the disk. Another interesting 
issue is the formation scenario of ULX in terms of binary 
star evolution. Supposing some ULXs arc highly accreting 
stellar-mass black holes, it is very natural to think that 
those ULXs are binary systems and the disk material is 
supplied by their companions. In the last part of this 
section, we comment on other plausible black hole mass 
estimation methods for ULXs. 
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4-1. Dependence on the Fitting Energy Range 

It is very important to note that the fitting results sen- 
sitively depend on the fitting energy range. In our present 
study, we fixed the fitting energy range to be 0.3 - 10 keV 
which is normally used in the study of BHBs as well as 
ULXs. Therefore, it is interesting to know how the fitting 
results will change if we choose a different fitting energy 
range. It will be necessary at some later time when the 
observations could provide us with data of a wider energy 
range. We extend our results to the higher energy part of 
the spectrum, so that we can also predict some features 
that had not been made available much by observations. 
On the other hand, we do not consider the lower energy 
part, lower than 0.3 keV, since there remains some diffi- 
culties in the analysis of this energy part due to the effect 
of the interstellar medium. 

We, therefore, performed the same procedure as what 
we explain in section 3 for a larger fitting energy range, 
0.3 - 30 keV. We summarize the fitting results in table 6. 
We obtained higher Xj n and lower r; n ; as a consequence 
we have larger mass correction factor, ranging from 1.4 - 
1.9 (m = 10), 1.4 - 2.0 (m = 30), and 1.4 - 2.1 (m = 100). 
To conclude our discussion in this subsection, we need 
a mass correction, m/m^ ~ 2 at most, for fitting with a 
wider energy range, which can be use as a reference when 
the observational data are available. 

4-2. Spectral Hardening 

The spectral hardening factor plays an important role 
in estimations of the black-hole mass, particularly when 
it is based on the X-ray spectral fitting. As mentioned 
in section 3, the spectral hardening factor is usually set 
to be 1.7 for the standard accretion disk. This value was 
obtained from a numerical calculation of the vertical struc- 
ture and the radiative transfer of an accretion disk around 
a Schwarzschild black hole by Shimura and Takahara 
(1995), in which a geometrically thin disk assumption is 
used. Their results imply that the X-ray spectrum of a 
standard disk is affected by Comptonization, and the lo- 
cal flux can be approximated by a diluted blackbody flux 
(see subsection 2.2), with k~ 1.7 (see also Czerny & Elvis 
1987; Ross et al. 1992). 

We are aware that the bound-free opacity may still play 
an important role, and that it is important to address this 
issue properly. In our present study, however, we did not 
include the effect of bound-free opacity, and we simply 
mimicked the deviation from the blackbody spectra by 
a diluted blackbody approximation without considering 
the details of the physical mechanism that may produce 
it. In fact, we also think that it is even possible that 
the spectrum at high-energy part may experience soften- 
ing as well (Ohsuga et al. 2003; Suleimanov & Poutanen 
2006). In order to address these issues properly, we need 
to work with radiative transfer, which is beyond the scope 
of our present study. Regarding the effect of the bound- 
free opacity on the spectral shape, it will be necessary for 
the readers to consult other papers (e.g., Suleimanov & 
Poutanen 2006). Meanwhile, we continue our discussion 
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Fig. 7. Luminosity vs. temperature diagram for various m (m = 10 in black '+', m = 30 in black 'x', and m = 100 in black '*') 
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are calculated from analytical solution of standard disk and red lines are calculated from standard disk relations (see text, section 
3.3). 



in the context of the diluted blackbody approximation in 
which spectral hardening factor plays a key role. 

Merloni et al. (2000) analyzed the spectral hardening 
factor used to correct the fitting results of the DBB model. 
They use a self-consistent model for the radiative trans- 
fer and the vertical temperature structure in a standard 
disk to simulate the observed disk spectra. In addition, 
they also take into account the gravitational redshift and 
transverse Doppler effects. They fitted this model to the 
DBB model, and found that the spectral hardening fac- 
tor is not constant; the spectral hardening factor becomes 
higher when the accretion rate and/or coronal activity is 
high. The range of the varying hardening factor is within 
1.7 < k < 3. Davis et al. (2006) also show significant spec- 
tral hardening when a disk becomes effectively optically 
thin. 

In terms of super-critical accretion model, Kawaguchi 
(2003) showed that spectral hardening by electron scatter- 
ing is quite a sensitive, since T os /r a b s c< a 1 / 8 . Therefore, 
an increase of a results in an increase of the color tem- 
perature, and thus we get more spectral hardening. Note 
that all the spectral calculations mentioned above assume 
hydrostatic disk layers. 

In the present study, we set the spectral hardening fac- 
tor to be, k = 1.7, the widely used value for the case of 



sub-critical (standard) accretion due to the uncertainty of 
the k in the case of super-critical accretion. The choice of 
other value of k would certainly affect the spectrum shape, 
for instance, choosing k > 1.7 will cause the spectrum to be 
harder in the high energy band and the slope in the lower 
energy band would also be changed, k ~ 3 was suggested 
for high luminosity disks (Watarai & Mineshige 2003b). 
Since the derived black-hole mass from the fitting depends 
on k as M oc re 2 , a mass underestimation/overestimation 
is possible due to the choice of re. 

To summarize, there remain uncertainties in the spec- 
tral hardening factor of super-critical accretion flow. This 
will, in turn, introduce another source of uncertainty in 
our present mass estimation. A more complete description 
of spectral properties of super-critical accretion is awaited. 

4-3. Outflow of Material 

Through a 2-dimensional radiation-hydrodynamic sim- 
ulation Ohsuga et al. (2005) found high-speed outflow 
(with 10-20 % of speed of light) from the center of super- 
critical flow with a rate of ~ Le/c 2 . This nearly rela- 
tivistic motion will affect the photon spectra, whose effect 
depends on the inclination angle of the system (Heinzeller 
et al. 2006). 

In fact, Heinzeller et al. (2006) found that the aver- 
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age photon energy due to a shift of the frequency of the 
escaping photons by the relativistic Doppler effect in a 
face-on case (small inclination angle) is of a factor of 1.18 
higher than that of the edge-on case (large inclination an- 
gle). They also reported that the main features of the 
thermal component of the spectrum of super-critical ac- 
cretion flow are consistent with the slim disk, while there 
remain unresolved features in the high-energy part. This 
kind of simulation is very important to address the effect 
of the outflow in the observed spectrum. However, an ad- 
equate calculation resolution has not been available at the 
moment. 

In the meantime, as a qualitative study, we compare 
two different cases. The first case is a disk with advection 
without outflow (our case), and the second case is the disk 
with outflow without advection, which was first consid- 
ered by Shakura and Sunyaev (1973), as also discussed in 
Poutanen et al. (2007). In addition to the fact that both 
cases result in similar temperature profiles, T G ff oc r" 1 / 2 
(flatter temperature distribution than that of the stan- 
dard disk), we found another interesting fact. 

Both cases have a characteristic radius. In the advec- 
tion case, it is called the trapping radius, r tr , inside of 
which the diffusion timcscalc is longer than the accre- 
tion timescale (and thus the photon trapping becomes 
important). In the outflow case, it is called the spher- 
ization radius, r sp , at which the disk luminosity is close to 
the Eddington luminosity (Lipunova 1999). Interestingly, 
those two radii are comparable or, to be more precise, the 
spherization radius is twice as small as the trapping ra- 
dius: r tr ~ mr g /2 while r sp ~ mr g /4, if we set H/r ~ 1 for 
the advection case. This means that the advection effect 
becomes effective even before the outflow needs to set in 
to remove any excess of the gravitational energy caused 
by high mass-accretion rate inflow of the disk materials. 

Although it is necessary to address this issue more prop- 
erly, especially when both advection and outflow occur si- 
multaneously (see also Poutanen et al. 2007), we did not 
include the effect of the outflow in our present study, while 
assuming that the advection works effectively. 

4-4- On Black-Hole Binary Evolution 

The work in the disk model is indeed synchronized by 
the work in the field of black-hole evolution in binary sys- 
tems. We point out some studies on black-hole binary 
evolution in this section. 

King ct al. (2001), for example, proposed a mild X- 
ray beaming model for ULXs. Their model suggest that 
ULXs may represent the phase of thermal-timescale mass 
transfer that occurs in high-mass X-ray binaries. This 
phase is short-lived, but is extremely common in the evo- 
lution of X-ray binaries. The short lifetimes of high-mass 
X-ray binaries could explain the association of ULXs with 
the star-formation region. Moreover, in the context of 
a binary system, their model does not require IMBH to 
explain the nature of ULXs. 

Rappaport et al. (2005) suggest that the luminosity 
of ULXs can be explained by a stellar-mass black hole in 
a binary with a moderate mass companion. They follow 



the evolution of the binary before the black hole is formed 
until it emits the X-ray by the accretion process onto the 
black hole. They found that luminosities of the order 
of ULX luminosities can be produced if the Eddington 
limit is allowed to be exceeded by about a factor of ten. 
Interestingly, this clearly supports super-critical accretion 
flow onto a black hole as a possible mechanism for ULXs. 

4-5. Other Black Hole Mass Estimation Methods 

Despite being very useful, some uncertainties remain 
concerning mass estimation through the X-ray spectral 
fitting as a result of its attachment on some assumptions. 
Therefore, mass estimation from other independent meth- 
ods are eagerly awaited to settle the dispute over the 
black-hole mass of the ULXs. It also remains possible that 
ULXs cannot be explained by a single theoretical model. 

Despite being widely used when the dynamical mass 
estimation method cannot be applied, the X-ray spectral 
fitting is not the only method to estimate the black-hole 
mass. Indeed, other independent mass estimation is de- 
manded to settle the dispute over the black hole mass of 
the ULXs. One of the other independent mass estima- 
tions is by double peak quasi-periodic oscillation (QPOs) 
in a 3:2 ratio to be discovered in the time variability of 
ULXs (Abramowicz ct al. 2004). These QPO frequen- 
cies scale inversely with the mass of the source, which was 
first shown for microquasars by McClintock and Rcmillard 
(2003), and more recently for low-mass Seyfert galaxies 
by Lachowicz et al. (2006). For ULXs, in general, it 
has not been detected yet. Recently, Strohmayer et al. 
(2007) reported possible double peak QPOs in NGC 5408 
X-l which happened to have 4:3 ratio. Unfortunately, the 
scaling law of the upper frequency QPO and the mass is 
not obvious for QPOs other than that with a 3:2 ratio. 

Other scaling laws are found in X-ray binaries and ac- 
tive galactic nuclei (AGNs). It was found that the ratio of 
the radio to X-ray emission depends strongly on the black 
hole mass (Falcke et al. 2004). If ULXs are in the same 
family with our Galactic X-ray binaries in the very high 
state, it is expected that they could also be radio tran- 
sients. In this regard, Kording et al. (2004) conducted a 
radio-monitoring study of ULXs in nearby galaxies with 
the VLA. Kording et al. (2005) found that generally ULXs 
do not emit steady-state radio emission above radio pow- 
ers of 1.5 x 10 17 W/Hz, and thus radio loud ULXs is still 
awaited to be found. 

Pooley and Rappaport (2005) suggest that X-ray and 
optical eclipses can be used to distinguish the two possi- 
ble black hole masses. ULXs with stellar-mass black holes 
should exhibit at least twice as many eclipses as would 
IMBH systems. They proposed X-ray monitoring cam- 
paigns for a few days to detect eclipses. In the optical 
range, monitoring ~ 10 ULXs twice a night each for a 
couple of months by using ground-based telescope should 
also work. 

Meanwhile, we can merely rely on X-ray spectral fitting 
which was found to be quite a reliable method in the case 
of Galactic black-hole binaries. 



No. ] 

Table 6. Same as table 4 but the energy range for fitting is 
0.3 - 30 keV. 



m 


m 


p 


kT~ in 




m/m x 


L/L E 


10 


32 


0.58 


1.39 


2.12 


1.39 


1.25 


10 


100 


0.53 


1.79 


1.66 


1.78 


2.10 


10 


320 


0.50 


1.96 


1.59 


1.86 


2.78 


10 


1000 


0.48 


2.04 


1.57 


1.88 


3.16 


30 


32 


0.57 


1.07 


2.08 


1.42 


1.26 


30 


100 


0.51 


1.41 


1.56 


1.90 


2.11 


30 


320 


0.48 


1.54 


1.49 


1.99 


2.79 


30 


1000 


0.47 


1.60 


1.48 


2.00 


3.17 


100 


32 


0.55 


0.80 


2.06 


1.40 


1.26 


100 


100 


0.49 


1.07 


1.48 


2.00 


2.11 


100 


320 


0.47 


1.17 


1.41 


2.09 


2.79 


100 


1000 


0.46 


1.21 


1.41 


2.09 


3.18 



5. Conclusion 

We have analyzed the assumption r; n = 3r g for black 
hole mass estimation for the case of super-critical accre- 
tion disk. 

We fitted the slim disk spectra for various parameter 
combinations (to, m, and i), with the extended disk black- 
body model. We fixed a = 0.01 in the disk structure cal- 
culation, and chose re = 1.7 for the spectrum calculation. 
To estimate the black hole mass, we adopt £ trans = 0.353. 

As expected from a theoretical point of view, we ob- 
tained p- values that deviate from 0.75, the value expected 
for the standard disk model. Moreover, we found that by 
adopting, the assumption r- m = 3r g for a slim disk, the 
estimated black hole mass should be corrected. The cor- 
rection factors are 1.2 - 1.3 and 1.3 - 1.6 for m = 10 and 
m = 30, respectively. Since the corrections are small, we 
can safely conclude that the black holes in ULXs which 
Vierdayanti et al. (2006) analyzed are stellar-mass black 
holes with mass < 100M Q . 
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